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Abstract 

We calculate the thermal flux of low-energy solar neutrinos and antineutri- 
nos of all flavors arising from a variety of neutrino pair processes: Compton 
production (including plasmon-pole diagrams), neutral current decay of ther- 
mally populated nuclear states, plasmon decay, and electron transitions from 
free to atomic bound states. The resulting flux density per flavor is significant 
(10 8 — 10 9 /cm 2 /sec/MeV) below ~ 5 keV, and the distributions fill much of 
the valley between the high-energy edge of the cosmic background neutrino 
spectrum and the low energy tails of the pp-chain electron neutrino and terres- 
trial electron antineutrino spectra. Thermal neutrinos carry information on 
the solar core temperature distribution and on heavy flavor neutrino masses 
for m u or m Vr > 1 keV. The detection of these neutrinos is a daunting but 
interesting challenge. 



A great deal of effort has been invested in detecting the solar electron neutrinos produced 
as a byproduct of solar fusion jy. This letter focuses on another source of solar neutrinos, 
those produced by charged- and neutral-current thermal processes operating in the solar 
core. While the mean energy of these neutrinos is low, on the order of the solar core 
temperature kT c ~ 1.3 keV, and their total flux quite modest, their flux density at earth is 
significant, ranging up to ~ 10 9 /cm 2 /sec/Me V per flavor. Thus, above cosmic microwave 
energies e v > 10~ 2 eV, this is the dominant source of low-energy heavy-flavor neutrinos, 
and of electron neutrinos below ~ 5 keV, where the low-energy tails of the solar fusion 
and terrestrial radioactivity neutrino spectra are encountered. In this letter we describe the 
processes producing such neutrinos and comment on the physics that might be extracted 
from their flux. 

The pair neutrino mechanisms we consider are illustrated in Fig. 1. The first is Compton 
production, a mechanism we find dominates the high-energy tail of the thermal distribution. 
The second is Compton production with a plasmon intermediate state, a process we believe 
has not been considered previously. The third is the decay of a plasmon into a vv pair. 
Both the ordinary Compton and plasmon processes are familiar neutrino cooling processes 
in stars 0. The fourth is the decay of a thermally populated excited nuclear state by pair 
emission, a process that can contribute at solar temperatures only in the case of an Ml 
transition at an anomalously low energy. Bahcall, Treiman, and Zee || identified the 14.4 
keV Mossbauer transition in 57 Fe as the principal contributor in the sun. The fifth process 
is pair emission as an electron makes the transition from a free state to a bound atomic 
orbital, the analog of an electromagnetic process known to be quite important to the solar 
opacity |§. 

The Compton process, as well as the other three leptonic mechanisms mentioned above, 
are governed by the weak Hamiltonian 



The couplings for heavy flavor neutrinos are gy = 1 — 4 sin 2 9 W and g& = 1 , while for electron 
neutrinos one obtains, after a Fierz transformation of the charged current interaction, the 
effective couplings gy = —(1 + 4 sin 2 6 W ) and g& = - 1- The Compton rate per unit volume 
is then given by 
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FIG. 1. Representative diagrams for the various thermal neutrino pair processes considered 
here: a) Compton process; b) plasmon pole contribution to the Compton process; c) transverse 
plasmon decay; d) nuclear Z° emission; and e) pair production in free-bound atomic transitions. 



The chemical potential //, which is determined by the electron density and temperature, is 
~ —2.1 keV at the solar core, so that the electron gas is slightly degenerate. 
The spin sum yields 

{m} e 

where the neutrino tensor is 

AT = 8[(pX + fift - w ■ PvgT) + i^ VvpVVc \. (5) 

After making a nonrelativistic reduction the electron tensor is found to be 

B w = J2 { 2 (9v + 9a){-^^ + <W • 5)(-we„ + 5 M e ■ q) + 2g 2 A [k fl k u - (e • qfg^]} ■ (6) 
H 

The remaining evaluation is simplified by the fact that T <^ m e , allowing one to treat the 
electrons in a static approximation. They absorb momentum but no energy, uj ~ u 7 = e u +e p , 
which leads to 
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where = 2/(27r) 3 / <i 3 pn e (p)(l — n e (p)) . 

The plasma produces an interesting modification of this result, illustrated in Fig. lb. The 
long-range Coulomb force between charges in the plasma produces a collective plasma oscil- 
lation, the longitudinal plasmon, characterised by the plasma frequency u p i. The transverse 
plasmons are photon modes that develop an effective mass equal to u p i (% = c = = 1). 

Neglecting the small contribution due to ionic charges, the plasma frequency is given by 

„ Aim e a . .,, i at / 

CfcCi = , where a is fine-structure constant, m e is the electron mass, and n e = NAP/Ue 

m e 

is the electron number density. Here is Avogadro's number, \x e = 2/(1 + X) is the 
electron's effective molecular weight, and X is the mass fraction of hydrogen. At the center 
of the sun, uj p \ ~ 0.28 keV. In the classical limit the dispersion relation for the transverse 
plasmon is 

u 2 =u 2 pl (l + q 2 /k 2 D ) + q 2 , (8) 

while that for the longitudinal plasmon is (0 < \q \ <oJ p i) 

u 2 = u 2 pl (l + 3q 2 /k 2 D ). (9) 

Here u; 7 is the energy of the plasmon and q is the momentum. Since the Debye wavenumber 

kr> = J e is much greater than uj p i, the plasmon energy is always very close to the 
plasma frequency. 

The plasmon contributions to the Compton process involve both the vector-vector po- 
larization n^] 7 and the vector-axialvector polarization n^ 4 . The terms involving UY^ are 
negligible, suppressed by (a^-) 2 in the rate. While the simple Compton and pole ampli- 
tudes interfere, this interference can be ignored since the pole contribution is so narrowly 
peaked at kinematics defined by q 2 ~ uj 2 v It follows that the pole contribution simply adds 
to Eq. (6) the following term 

B£ = 2<#[(e-V(e-% - 9ila g pv \ E("^ CT + 5 a0 e ■ q){-ue» + 5 p0 e ■ q), (10) 

where (e^ 1 )^ = g^ v + D pa Il a u is the inverse of the usual dielectric tensor of the plasma, 
D^u being the full photon propagator. (The zeros of the determinant of the dielectric tensor 
determine the plasmon modes.) As this contribution is proportional to g 2 ,, it is numerically 
significant only for electron neutrinos. 

The net contribution of the plasmon pole can be evaluated by integrating over a sharply- 
peaked off-shell region determined by the small width of the plasmon. The plasmon width 
is zero in the random phase approximation (decay into one-particle one-hole electron states 
is forbidden for timelike q 2 ). Thus one must go to the next order: we take the width from 
the two-particle two-hole calculation of DuBois et al [§ . As we are primarily interested in 
neutrino energies above u p i, the longitudinal mode can be ignored. A somewhat tedious 
calculation then yields the additional contribution to the Compton rate from transverse 
plasmons 
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dt " A%^mf v J " " 7 ^Ime t 
The imaginary part of the transverse dielectric function is 



lme t = Y^£^[^K)/3 + (12) 
where the plasma parameter A = k^/n e , and 

I(u; 7 ) ~ ^ ^ exp[-x - ^/(16T 2 x)]. (13) 

These expressions are valid for cu 7 3> and \q\ <C ~ 5.5 keV. The lower bound on cj 7 is 
important: plasmon contributions are large at low energies, so that later numerical results 
should be viewed with caution in this region. The constraint for large \q\ is not a concern. At 
large |g| and thus large e u + e p , the pole contribution is small and thus uninteresting relative 
to other processes: the pole diagrams become 0(a) corrections to the vector-coupling part 
of the simple Compton process. 

We note that the analogous Compton plasmon pole calculation for very hot or dense 
stars would require a finite-temperature relativistic calculation for the imaginary part of 
the polarization insertions beyond the random phase approximation. We believe no such 
calculation has yet been performed. 

The contribution of direct transverse plasmon decays into vv is illustrated in Fig. lc. 
Note that this process and the Compton plasmon pole process are entirely distinct: the 
plasmon is off shell in the pole process. The decay rate per unit volume is given by 

d -w = /" w w *£)■ W<<M( » +p ° ~ q) g w • (14) 

where J2{e} denotes the summation over the two polarizations and 

\M 1 \ 2 = A^(T, a e a )(Y v ^T , (15) 

with T^ a the effective photon-neutrino coupling || and e a the photon polarization four- 
vector. The neutrino tensor A^ u is given by Eq. (5), while can be determined from the 
VV and VA polarization insertions. The latter is again negligibly small. We find 

r^e a = -^L=g v Il t ( q )(0,er, (16) 

where the transverse VV polarization insertion can be approximated by II t (g) ~ uj 2 v The 
final result then reduces to 

dN G 2 g 2 v 6 f 1 1 ( 2 , 2W , , 17 v 



dt 128vr 4 a pl J u 2 e^' kT 
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where u 1 = e v + e p and the integrations extend over all neutrino energies satisfying cj 7 > uj p g. 
We have ignored the contribution of longitudinal plasmons, which again contribute only if 

£u + ~ LO p i. 

Neutral current decays of excited nuclear states can produce pair neutrinos (Fig. Id). 
The importance of a given element depends on the probability that the excited state will 
be thermally populated, on the element's solar mass fraction, and on the strength of the 
isovector axial (Ml) transition matrix element. As noted in Ref. ||, the important metal 
is 57 Fe because of the low-lying 3/2~ — > l/2~ Mossbauer transition, despite the low solar 
abundance of this isotope (mass fraction of 3.26 -10~ 5 [0]). We can express || the pair 
neutrino rate per flavor in terms of the known gamma decay rate uj 1 = 7.41 • 10 5 /sec (thereby 
eliminating most of the nuclear physics uncertainty) 



dt 1 + 5 2 b< 7T 3 a W$ \ji 1 + (ii -l)i3-r ] ) 2 Jo 

where is the number density of 57 Fe nuclei in the solar core, Wq = 14.4 keV is the 
3/2~ — > l/2~ transition energy, 5 = 0.0022 is the E2/M1 mixing ratio, /Iq = 0.88 and /ix = 
4.71 are the isoscalar and isovector magnetic moments, and Fa = 1.26 is the axial vector 
coupling constant. The nuclear-structure-dependent terms 

<Jf\\T£j(i)T 3 (i)U> ,, ogo 

< J f \\Y,l 1 f(i)T 3 (i)\\J i > 

3= <jf\\Ttim\Ji> 119 (19) 

are taken from the shell model calculations of Ref. ||. Note that the maximum rate for 
Eq. (18) occurs when Wo = 5T, or about 6.5 keV in the solar core, so the 57 Fe transition 
would be more effective in somewhat hotter stars. 

The final neutrino pair process we consider is Z° emission in the transition of a continuum 
electron to a bound atomic orbital (Fig. le). Although such free-bound transitions to atomic 
states in Fe and other metals are known to be an important contributor to the solar opacity 
0, relatively little attention has been paid to the analogous weak process. Early work was 
flawed by the use of plane waves for the continuum state: the resulting nonorthogonality 
of the initial and final wave functions leads to an overestimation of the rate by a factor 
proportional to m e /T. More recently the process was treated correctly by Kohyama, Itoh, 
Obama, and Mutoh 0, who tabulated total cooling rates for a variety of metals using the 
dipole approximation (which we also adopt) and the transition to the n = 1 i = atomic 
bound state. 

The rate for capturing into Bohr orbits described by (n, £) is 
dN r 2 

- = N(Z) fin, / d y_____ Vp a ip) (20) 

where e p = e v + e p — E b (n,£) is the continuum electron's energy, E b (n,£) is the (positive) 
Bohr orbit binding energy, N(Z) is the number density of the atoms in the solar core, v p 
is the electron velocity, a(p) is the cross section for capture of an electron of momentum p 
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on the ion (averaged over initial electron spin and summed over all final lepton spins), and 
f(n,£) is a blocking factor correcting for the incomplete ionization of the atomic shells. (In 
the solar interior, the only atomic shells not completely ionized are those with low n and 
high Z.) 

The differential cross section can be readily evaluated in the dipole approximation, as 
discussed in Ref. |J. The result is 

ff (P) = -l-U 7 ft<" *\r\p £-l)\ 2 + (£ + l)\(n £\r\p £ + 1)| 2 ) 
x f 1 dx{\ - xfx 2 \gl{2x 2 - 2x + 3) + 2g 2 A (5x 2 - 5x + 3) 

(21) 



±— 0^a(16z 3 - 24x 2 + 14x - 3) + (— ) 2 c£(28x 4 - 56x 3 + 46x 2 - 18x + 3) 



where q = e u + E b (n,£) = e u + e p , x = e u /q , and where the radial integrals are defined 
with unit normalization for the incoming continuum electron wave function of momentum p. 
The ± sign on the third term distinguishes the v (+) and v (-) spectral distributions. The 
radial integrals between hydrogenic bound states and spherical components of the incoming 



Coulomb wave of momentum p can be evaluated in closed form jlOfl . The integration over 
the neutrino spectrum can be readily performed to give 

<*> - + & + m - ^ 

where M 2 is the matrix element of Eq. (21). The first two terms agree with the corresponding 
calculation of Ref. 0. In using Eq. (21) to evaluate the free-bound neutrino spectrum, we 
ignore the terms proportional to qo/m e , which are clearly quite small at solar temperatures. 
This eliminates the VA interference term that distinguishes the v and v distributions. 

The qualitative features of free-bound transitions can be understood readily. As the 
four-momentum transfer to the neutrinos is timelike, the corresponding momentum transfer 
|<f| < <?o ~ %T ~ 2 keV, so that \q\a /Z ~ 0.5/Z, where a /Z is the Bohr radius of the 
element of charge Z. Thus the dipole approximation used above is somewhat marginal for 
hydrogen, but quite reasonable for heavier species. In our calculations we took abundances 
for H, He, C, N, and O from standard solar model results, and adopted Cameron abundances 
for heavy metals (Ne, Mg, Si, S, Fe). The results (see Fig. 2) show that hydrogen accounts for 
about 1% of the captures, despite its high abundance: free-bound transition are dominated 
by the heavier metals. Thus the dipole approximation is very well justified. 

The importance of the heavier metals results from the fact that the electron's appreciable 
momentum, p e ~ y / 2m e (|T) ~ 45 keV, must be absorbed by the bound state, so that 
compact high Z atomic states containing such high momentum components are favored. 
The dipole matrix elements depend on the parameter 77 = Z/a p e ~ 0.08Z. For small rj 
(thus small Z), the transition probabilities vary as r\ 21 and strongly favor capture into s 
states. As long as 77 is small, the n — 1 state dominates the capture, accounting in hydrogen 
for 97% of the total at the peak of the free-bound neutrino spectrum. However this is 
not the case with high Z ions, when 77 becomes comparable to n: the n — 1 transition in 
Fe accounts for only ~ 60% of the total cross section. In our calculations we summed all 
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FIG. 2. The total free-bound v e flux density (solid line) and the separate contributions of Fe, 
Ne+Mg+Si+S, the CNO elements, He, and H are shown as a function of the neutrino energy. 
Ignoring corrections of order qo/m e , the v e results are identical. The heavy- flavor neutrino flux 
density is about a factor of 0.46 smaller. 

transitions through n = 5, which we found necessary for achieving 0.1% accuracy in the case 
of Fe. 

Figure 2 shows that the heavy a-stable metals, Ne, Mg, Si, and S, dominate the free- 
bound rate near the peak of the distribution, but give way to Fe at large e v . This reflects 
the kinematic advantage of Fe because of its larger atomic binding energy. At very low 
energies the CNO elements are the largest contributor. One might hope that the thermal 
solar neutrino flux could be exploited to probe the metallicity of the solar core, testing 
one of the assumptions of the standard solar model. As the sun is believed to have been 
highly convective at the onset of the main sequence, and thus homogenized, the standard 
solar model equates the initial core metal abundances to present-day surface abundances. It 
would be nice to verify this assumption experimentally by detecting the free-bound neutrino 
flux. Unfortunately we will see below that the free-bound flux is generally hidden by other, 
stronger thermal neutrino processes. 

The contributions to the terrestrial neutrino flux from solar thermal neutrino processes 
were obtained by summing the contributions from the five processes described above and 
depicted in Fig. 1. The rates were folded with the standard solar model density and temper- 
ature profiles which were also employed in calculating the electron chemical potential 



as a function of radius. Integrals were done over the solar core with a sufficient number of 
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Gaussian mesh points to achieve a numerical accuracy of at least 0.1%. The distributions of 
heavy metals not taken from the standard solar model (Ne, Mg, Si, S, Fe) were assumed to 
be uniform throughout the sun: no attempt was made to estimate effects due to diffusion. 

Figure 3 gives the resulting v e jv e and heavy-flavor flux densities. (Note that certain 
thermal processes, such as bound-bound transitions, will contribute but have not been es- 
timated: as bound-bound transitions are known to be relatively unimportant to the solar 
opacity, we assumed they would also be small in the weak case. This process can be con- 
sidered an omitted correction due to incomplete atomic ionization in the solar core.) The 
figure shows that the simple Compton process dominates the high-energy tails of both the 
v e j v e and heavy-flavor spectra, e v > 5 keV. In the heavy neutrino spectrum, the free-bound 
process dominates below 2 keV. Clearly the only "window" on core metallicity is then a 
tiny one, involving the enormously difficult task of observing very low energy heavy-flavor 
neutrinos. In the case of electron neutrinos, the pole contribution to the Compton process 
dominates below 5 keV. (The pole Compton and plasmon rates are proportional to g v and 
thus are very much smaller for heavy flavor neutrinos.) One concludes that the entire v e 
spectrum - this is the species most likely to be detected some day - is governed by purely 
leptonic thermal processes that can be readily calculated from the core temperature distri- 
bution and electron density. Thus it is information on the core temperature that is encoded 
in this flux. 

In Fig. 4 we superimpose the spectrum of thermal neutrinos, generated by summing all of 
the processes discussed here, on a familiar graph [12| of sources contributing to the neutrino 
flux density at earth. The thermal contributions are significant. The peak of the thermal 
solar v e (or z/ e ) flux density occurs at low energies, where the dominate production mecha- 
nisms are plasmon decay and the plasmon pole contribution to the Compton process. The 
peak flux density is ~ 10 9 /cm 2 /sec/Me V, a value somewhat less than 1% that achieved at 
the peak of the pp solar fusion neutrino distribution and about an order of magnitude larger 
than the terrestrial v e peak. The heavy-flavor flux density reaches its maximum at ~ 2 keV, 
achieving a peak value about an order of magnitude below that found for electron neutrinos. 
It is the only important heavy-flavor contribution in Fig. 4 above cosmic background ener- 
gies. The thermal v e j v e spectrum is the dominate feature above cosmic background energies 
and below ~ 5 keV, where the solar fusion v e and terrestrial v e distributions take over. The 
integrated v e and fluxes are 2.1 • 10 6 /cm 2 /sec and 3.0 • 10 5 /cm 2 /sec, respectively. 

In using Fig. 4 it is important to remember that in our evaluation of the Compton pole 
process, and in our neglect of longitudinal plasmon decay, we have made approximations 
that should break down near and below t v ~ uj v \. Thus these contributions to the low-energy 
results (< 0.3 keV) in Fig. 4 are a naive extrapolation, and should be viewed with caution. 

The thermal neutrinos are unimportant to solar evolution as they account for a neglible 
fraction of the energy loss (~ 0.001%). However, the possibility that they could be exploited 
to check the solar core temperature is intriguing. The rate of the Compton process varies 
as T 7 . Thus a 5% change in the core temperature - a variation that would significantly 
reduce the 8 B solar neutrino flux - would be expected to change the thermal neutrino flux 
by about 30%. The high-energy tail of the distribution would be considerably more sensitive 
to temperature variations, of course. This cross check on the standard solar model might 
be interesting because it is independent of nuclear cross sections, on composition gradients, 
or any other detailed aspects of the model. Another nice property of such a check is that 
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Neutrino Energy (keV) 

FIG. 3. Neutrino flux densities from free-bound transitions (long dashes), the Compton pair 
process (intermediate dashes), the pole contribution to the Compton pair process (short dashes), 
57 Fe nuclear decay (dashed-dot), and transverse plasmon decay (dotted). Totals (solid line) are for 
a single flavor. 



the thermal v e flux is left unaffected by the simplest MSW neutrino oscillation scenario [13 



z/ e s do not experience a level crossing because the effective mass of the u e has the wrong 
sign, assuming the usual mass hierarchy where the electron type is the lightest flavor. This 
contrasts with the case of solar fusion u e s: in principle, the comparison of the pp, 7 Be, and 
8 B neutrino fluxes is very sensitive to the core temperature, but in practice the fluxes appear 
to reflect neutrino oscillations, making the "thermometer" somewhat harder to read. Thus 
a second thermometer with quite different sensitivity to new physics (both in the MSW 
mechanism and in Sm 2 /e u ) could be quite valuable. 

In Fig. 5 we show where u e s of different energy are produced within the sun. Clearly the 
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FIG. 4. Natural neutrino sources. The terrestrial v e flux and continuous flux of extragalactic 
supernova neutrinos of all flavors are from Krauss et al. 12]. The solar (fusion) v e flux is the 



standard solar result of Bahcall et al. [11]. The thermal solar neutrinos are for a single flavor. 



spectrum contains a great deal of information on the temperature distribution within the 
sun. 

These remarks are made because the most likely opportunity for measuring the thermal 
neutrino spectrum is a process that depends on flux density, not on total flux, and which 
samples that flux at a precise energy, the resonant reaction 

u e + e- + (A,Z)^(A,Z-l). (23) 

This reaction has been discussed previously in connection with terrestrial v e sources ||12|| . 
Cross sections can be large in high Z atoms, where the electron overlap with the nucleus 
is favorable. Because nuclear level widths are very narrow, this process samples the v e flux 
density at a discrete energy. The are several possible candidate transitions with energies 
between 2 and 20 keV. (One that has been studied in connection with neutrino mass mea- 
surements is the decay of long-lived 163 Ho to 163 Dy, which has a positive q- value of less than 
3 keV: either a neutrino mass or v e inducement of electron capture alters the atomic orbits 
that participate in the capture.) 

The heavy-flavor neutrino flux also contains interesting information: if the existence of 
this flux were established, it would immediately impose kinematic mass limits of ~ 1 keV on 
the and v T . Unfortunately there is no obvious possibility for measuring these species. The 
problem could well prove as difficult as in the case of the cosmic microwave neutrinos, where 



existing experimental bounds exceed the expected flux by about 15 orders of magnitude [14 
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15.0 




FIG. 5. The radial distribution of neutrino pair production within the sun for three selected 
energies, 0.33, 4.3, and 12.6 keV. Each distribution, integrated over R, is normalized to unity. 
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